Double- Weyl Topological Metals Stabilized by Point Group Symmetry 
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We study the properties of 2D insulators and 3D metals which possess both n-fold rotational 
symmetry and broken time-reversal symmetry. For 2D insulators, we show how to determine the 
Chern number up to a multiple of n by evaluating the eigenvalues of the rotational operator at high 
symmetry points within the Brillouin zone. For 3D metals these eigenvalues also determine the 
existence, number and type of the gapless band crossings. We apply this theory to show that [001]- 
axis ferromagnetic HgCr2Se4 in 3D possesses both a Ca symmetry-protected quadratic (in certain 
directions) band crossing (dubbed double- Weyl) as well as a mirror-symmetry protected nodal line 
separating regions of different electrical polarization which give rise to a multiferroic behavior. If the 
direction of the ferromagnetism is shifted away from the [001]- to the [lll]-axis, a combined Cq and 
mirror symmetry reveals a difference of two in the Chern number of different crystalographic planes, 
giving rise again to a topological metal. The possible quadratic-crossing is, however, not protected 
by any symmetry, and small spin-orbit coupling terms break it into four linear Weyl fermions. 

PACS numbers: 



The discovery and classification of topological states 
of matter beyond those found in time reversal invariant 
(TRI) materials [ll-lsl is vital to understanding the full 
spectrum of possible topological insulators. In principle, 
it is possible to find many more topological systems with 
distinct physical observables characterized by the sym- 
metries they preserve. Examples of topological systems 
possessing symmetry beyond TRI and charge conjuga- 
tion are systems invariant under point group symmetries 
(PGS)|9l4ll|. Their topological classification does not in- 
volve TRI, and, therefore, these systems may encompass 
a larger class of materials including magnetically ordered 
ones useful for potential applications [Til [l3| . 

Most recently, 3D topological semimetals have been 
proposed to exist Jn a variety of materials such as 



cold atom systems [1J|, multilayer topological insulator 
systems [TBI, flBj and pyrochlore iridates[17, 18]. Unlike 
in a topological insulator, the band structure of a topo- 
logical semimetal exhibits gapless points in the Brillouin 
zone (BZ), in the close proximity of which the effective 
Hamiltonian is that of a 3D Weyl fermions with linear 
dispersion. In a translational invariant system, these 
gapless points are protected from opening a gap through 
infinitesimal transformations of the Hamiltonian; these 
points act as monopoles (vortices) of 3D Berry curva- 
ture, any closed 2D surface surrounding them exhibits a 
unit Chern number, and they can be gapped only by an- 
nihilation with other Weyl points of opposite monopole 
charge. They can be present both in the presence of TR 
(the metallic state between a trivial and an inversion non- 
symmetric topological insulator with TR is a topological 
metal with 4 Weyl points) and in its absence (pyrochlore 
irridates). Without any other symmetries, Weyl points 
are the generic topological band-crossings in 3D. 



In this letter, we show theoretically that a series of new 
3D topological metals can exist when PGS are present 
in a material. The metals exhibit double (quadratic in 
certain directions) or multiple-crossings at gapless points 
in the BZ. These crossings, which we dub "double Weyl 
points" , are protected from splitting into Weyl fermions 
by the PGS. We first investigate 2D Cn rotation sym- 
metrical topological insulators with broken time-reversal 
symmetry, and find that their Chern numbers can be 
determined (mod n) from the Cn eigenvalues of the oc- 
cupied bands at high- symmetry points in the BZ. We 
then show that a 3D topological metal with PSG whose 
Hamiltonian, a function of the third momentum, interpo- 
lates between two 2D topological insulators with differ- 
ent Chern numbers can exhibit double gapless points of 
quadratic (or higher order) touching. We analyze the re- 
cently proposed semimetal HgCr2Se4[l^ and find it is a 
double- Weyl topological semimetal stabilized by C4 sym- 
metry. As a result of changing mirror symmetry eigen- 
values, we find that HgCr2Se4 also exhibits a topological 
nodal line separating regions of different charge polariza- 
tion. More generally, we find that the existence, number, 
and nature of the crossings in topological PGS materials 
may be determined from the eigenvalues of the symmetry 
operator at high symmetry points within the BZ. 

Consider a generic 2D system with n-fold rotation 
symmetry, Cn (lattice periodicity dictates that n = 
1,2,3,4,6), and translational symmetry. A simple ar- 
gument shows the existence of a relation between the 
Chern number (mod n) and the eigenvalues of the occu- 
pied bands under Cn symmetry at high symmetry points. 
As in the C2 case of [9|, llO| , starting from an atomic limit 
insulator and changing the hopping parameters, the sys- 
tem could undergo gap-closing and reopening transitions 



2 



(in 2D the two-band crossing co-dimension is 1) towards 
a topologically non-trivial insulator with nonzero Chern 
number. If a band crossing occurs at a generic point in 
the BZ, Cn rotation symmetry fixes the location of the 
other necessary n — 1 crossings, giving a change in Chern 
number of n (= mod n), but chang ing no Cji eigen- 
values. However, less than n band crossings with total 
Chern number change less than n are possible if they 
take place at high-symmetry points, but these are now 
accompanied by changes in the Cn eigenvalues of the oc- 
cupied bands. This argument (alternatively, one can use 
the monodromy method can be placed on analytical 
footing to obtain 



i-ifi^ = n ur)ux)UM), (1) 

n^occ. 

(-i)^exp(iyC) = n On{r)en{K)en{K'), (2) 

n^occ. 

(-l)^exp(i|c) = n Vn{r)0n{K)UM), (3) 

nGocc. 

where ^ V eigenvalues of C2, C3, C4, Cq 

rotation operators respectively, and F = (F = 1) for 
spinless (spinful) fermions. Alternatively, the symmetry 
eigenvalues of the conduction and valence bands fix the 
form of the effective Hamiltonian around a possible two- 
band crossing at a high-symmetry point. 

We begin by examining systems with C4 symmetry. If 
through a phase transition the Chern number changes by 
a multiple of four then the requisite band crossing may 
occur anywhere within the BZ. However, if the change is 
not a multiple of four, there must be one crossing at any 
one of the four high symmetry points in the BZ namely: 
r, X, Y and M. A crossing at X (or F), necessitates 
another crossing at Y (or X), as these two points are 
equivalent under C4. At the X and F, only the inversion 
operator C2 = Cl commutes with the Hamiltonian so we 
revert to the case with inversion symmetry of [ol, [lo| . If 
the two bands that cross at X or F have different in- 
version eigenvalues, the dispersion must be linear around 
X or y and the total change in Chern number must be 
with each point contributing a change of one. 
The situation is different at the stand-alone points V and 
M. We may discuss them separately, noting that each 
crossing contributes additively to the Chern number. At 
each point, there are two types of band-crossings between 
two bands with different eigenvalues: Case (i), one band 
has the C4 eigenvalue e*^^/^ or — e*^^/^, and the other 
g-2F7r/4 _g-zF7r/4 Cd^^e (ii), ouc band has the C4 



or one band 



eigenvalue e*^^/^ while the other — e*^^/^, 
has e~*^^/^ while the other — e~*^^/^. 

For Case (i), since C2 = Cl, we have a crossing be- 
tween two bands with opposite inversion eigenvalues. At 
the crossing point, the Chern number changes by one, 
but with only inversion it is impossible to tell a priori 
whether the change is positive or negative as P is only 




FIG. 1: (a), (b), (c): the three crossing scenarios that change 
the Chern number C{kz) of the 2D Hamiltonian H{kx, ky;kz) 
by two as kz increases from kz = to kz = 7v/2. In (d) the 
subsystem consisting of all the k's inside the cylinder has a 
non-zero electric polarization, with each {kx, ky) contributing 
e/2. The dotted red line in (d) is the nodal line topologically 
protected by x^-plane mirror symmetry. 



defined up to a gauge. C4 symmetry fixes the gauge 
and hence allows for further classification of whether the 
Chern number increases or decreases by one. When con- 
sidering Case (ii), the two bands have the same inversion 
eigenvalue and thus, in the basis given by the two cross- 
ing states, the C4 operator must be cr^e*^^/^ or ia^e*^^/^, 
corresponding to the case ^1 = = — e*^^/^ and 
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F7r/4 respectively. From here we 



take the former case C4 = e^^^/^a^ and the solution pro- 
gresses in a similar fashion to C4 = le^^^l^Gz- 

The most general form of an effective two-band Hamil- 
tonian is i^e//(k) = e(k) + (ii(k)cri.C4 (= Oz) rotational 
invariance C^^Heff{kxjky)C4 = Heff{ky^—kx) requires 
di,2{kx^ky) = —di^2{ky^—kx)^ which means di(k)^d2(k) 
cannot be constants or linear functions of k in the prox- 
imity of the r or M points. There are three pos- 
sible quadratic combinations of kx and ky: k^ -\- ky, 
kl — ky and kxky, with only the d-wave form allowed 
in (ii(k), (i2(k). Any linear combination of k'^ — ky and 
kxky is allowed in (ii^2(k). C4 invariance also implies 
dsikx^ky) = ds{ky, —kx), satisfied by a constant mass 
term m or by the s-wave k^ + ky (which we neglect for 
small k in comparison to m), resulting in an effective 
Hamiltonian of: 



eff 



(4) 



where k± = kx ^ iky, (j± = ± icjy and f,g are in 
general complex numbers depending on details of the 
model. A direct evaluation of the Chern number yields 
sign[m(|/| — 1^1)] for the upper band and — sign[m(|/| — 
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1^1)] for the lower band. A band crossing where m re- 
verses sign changes the Chern number of the occupied 
band by ±2 while at the same time switching the C4 
eigenvalue (= — sign(m)) of the occupied band at F or M 
from 1 to —1 or viceversa. 

We now bridge from 2D insulators to 3D semimetals 
and show how our PSG eigenvalue equations [3] can help 
us understand and predict the physics of realistic sys- 
tems. A 3D system H{kx,ky^kz) with an n-fold rota- 
tion axis (defined as the z-axis unless otherwise noted), 
can be represented as a /c^-dependent sequence of 2D sys- 
tems H{kx^ky]kz) with Cn symmetry. The C4 eigenvalue 
changing crossings discussed above at F, X, F, and M 
can now happen on high symmetry lines along z-axis that 
pass through these 2D high-symmetry points. With this 
understanding, we examine HgCr2Se4 using an 8-by-8 
tight-binding generalization of the k • p Hamiltonian ma- 
trix [l9| which exhibits 7r/2 (spin and coordinate) rota- 
tional symmetry about c-axis or C4 = ex.p{i{Lz-\-Sz)7r/4:): 

C4 H(^kx^ ky^ kz)C4 = H(^—ky^ kx^ kz)' (5) 

In addition, the system is also invariant under 3D in- 
version, P. As Mz = C4P, the system is also invari- 
ant under a mirror symmetry about the x?/-plane. If 
at a fixed /c^, the energy spectrum of H{kx^ky]kz) is 
gapped for all kx and ky, the Chern number on that 
plane, C{kz) may be conveniently calculated by applying 
_^C(M = Yl^^^^^^,{0,0,kzM7r,n,kz)Un,0,kz). The 
face-centered cubic lattice of HgCr2Se4 is invariant under 
a translation of niai +n2a2 +n3a3, where ai = (1, 1, 0)a, 
a2 = (1,0, l)a and as = (0, 1, l)a (a = 5.375A). A 
Bloch Hamiltonian has the reciprocal vector (in units 
of a~^)translational symmetry H{kx,ky,kz) = H{kx ± 
TT^ky ± TT, /c^ ± tt). As the Chern number in a 2D lat- 
tice Hamiltonian is invariant under any in-plane transla- 
tion in /c-space, we have Ci{kz) = Ci{kz + tt). There- 
fore if Ci{kz = 0±) ^ Cilkz = 7r/2), {0± = ± S 
with 6 being an infinitesimally small positive number, 
as the kz = plane exhibits a nodal line) then from 
kz = 0+ to kz = TT the Chern number must change twice 
at kz = kc and kz = 7T — kc {0 < kc < 7t/2). For HgCr2Se4 
in the presence of ferromagnetism along the [001] direc- 
tion, we compute the eigenvalues of C4 at high symmetry 
points on both kz = 0± and kz = 7r/2 and conclude 
that Ci{kz =0) = 4n + 2 and Ci{kz = 7r/2) = 4n. 
A direct calculation of the Chern number on a nearest 
neighbor hopping tight-binding model on the FCC lat- 
tice for HgCr2Se4 confirms our eigenvalue formula, giving 
Ci{kz = 0±) = -2, Ci{kz = 7r/2) = 0. While the Chern 
number change was shown by a first principle calculation 
in [l9|, the relation to the C4 symmetry is new. 

As C{kz) is not uniform, HgCr2Se4 is necessarily a 
metal. The most dramatic effect of the PSG is mani- 
fested in the nature of the crossing point (s). The crossing 
must happen along either TZ or XP in the first BZ of 
the body-centered tetragonal. If it happened anywhere 



else, by (^4 symmetry, three other points in the BZ would 
also experience the same type of band crossing, leading 
to a Chern number change (as a function of kz) which 
is a multiple of four. In principle, there can be three 
crossing scenarios for a Chern number change of two: (1) 
two linear crossings on FZ, as seen in Fig. [D^a), (2) one 
double- Weyl crossing on FZ, as seen in Fig. [IJb), or 
(3) one linear crossing on XP , as seen in Fig. [D^c). 
Along FZ, C4 commutes with 7^(0,0,/^^) so all bands 
can be labeled by their C4 eigenvalues. Therefore, we 
only need the eigenvalues of the two crossing bands to 
classify the type of crossing. The HgCr2Se4 bandstruc- 
ture, seen in Fig. Efa), tells us that, along TZ at some 
kz = kc^ two bands with ±e*^/^ C4 eigenvalues respec- 
tively cross each other (see inset of Fig. [3j) resulting in a 
double- Weyl point (as a function of kx^ky the Hamilto- 
nian takes the m = form of Eq[4]), seen in Fig. [Ifb). 
This is the main difference between HgCr2Se4 and other 
proposed systems thought to support linear dispersion 
nodes (Weyl nodes) |16l4l8| which are spatially segregated 
and can only accidentally come together to either anni- 
hilate or double- up. In HgCr2Se4, the double Weyl point 
can be thought of as two single Weyl fermions of same 
monopole charge stuck together and prevented from be- 
ing separated by C4 = e^^/^a^ symmetry. The double 
Weyl point is indeed stable to infinitesimal changes of 
the Hamiltonian. If a spatial physical cut is made par- 
allel to the x-direction, we expect the appearance of two 
Fermi arcs (see Fig. [3j) on the cross section, that origi- 
nate from one of the bulk double- Weyl points and end in 
the other. 

Another PCS available in many materials is mir- 
ror symmetry. The existence of a xy-imiioi plane 
implies, for a 3D system, that H{kx^ky^—kz) = 
MzH{kx^ky^kz)Mz' Hence, at every (kx^ky), the 
Hamiltonian H(kx^ky^kz) possesses ID kz inversion sym- 
metry which quantizes the electric polarization to be 
Pzikx^ ky) = or Pzikx^ ky) = e/2 depending on the for- 
mula (-1)2^^ = n^Gocc.C(0)C(^)Iio[. Changing values 
of Pz from to 1/2 implies going through a bulk phase 
transition. In HgCr2Se4 we find that at a point very 
close to kx = ky = 0, Pz{d^5) = 1/2 {dX kx = ky = 0, 
the system is gapless at kz = i/^c, hence we compute the 
polarization at kx^ky = 0^), while at another corner of 
the BZ, Pz{7r,0) = 0. If we consider any path connect- 
ing {S,6) and (7r,0), there must be a point on the path 
where the gap closes. As the path is arbitrary, the re- 
sult is a closed loop of nodes, or a line node in the BZ, 
due to the mirror symmetry. The nodal line can only ap- 
pear on kz = (the case here) or the kz = tt plane. For 
{kx^ ky) outside the nodal line, the electric polarization is 
quantized to zero and within the nodal line, every point 
contributes eL/2 along z-axis (as seen in Fig . [Hd)), 
where L is the length of the system along z. The bulk 
polarization becomes Pz = ^ 0.64/iC/cm^ where 

S'at is the area encircled by the nodal line on the kz = 
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FIG. 2: Dispersion and pseudospin configurations of the var- 
ious protected crossings between the conduction and valence 
bands in FM HgCr2Se4 with FM//[001] (a,b) or FM//[111] (c- 
f). (a) C4-protected double- Weyl node with kz — 0.43 in BZ; 

(b) Pseudospin configurations along equal energy contours of 
the valence band in (a) at ^ = -0.02, 0.10, -0.13eV (blue, 
green and red arrows); (c) and (e) Dispersions around the 
one and three (related to each other by C3 symmetry) Weyl 
nodes at kz — 0.43 and kz — 0.15; (d) and (f) Pseudospin con- 
figurations along equal energy contours in the valence band of 

(c) and (e) at Ei = 0.0317eV and E2 = 0.02eV respectively. 
The winding in (d) is opposite to that in (f), indicating the 
opposite monopole charge of the k'^ = k'y = Weyl node from 
that of the three Weyl nodes away from the VZ' axis. 



plane. While ID polarization is gauge variant, should one 
wish to use a different gauge, the same ID, /c^-dependent 
gauge would need to be applied to the entire kx^ ky pa- 
rameters. We note that the electric polarization in this 
system comes from the particular band structure topol- 
ogy of the ferromagnetic state, and thus, HgCr2Se4 may 
be referred to as a "topological multiferroic" . 

Eq. ([3]) also predicts a topological metal if the fer- 
romagnetism in HgCr2Se is now aligned with the [111]- 
direction (which we now denote as the new direction, 
in a system {x'^y'^z')). The existence of a topological 
metal in this configuration was first noted in [19| from 
LDA studies but the link to the PSG was not made. In 
this system, the PCS is Cq^Mz'^ i.e., a combined symme- 
try of a six fold rotation symmetry about z' and a mirror 
reflection about the x'?/'-plane. The translational sym- 



metry of the BCC reciprocal lattice determines the period 
oi C{k'^). In the new frame, the invar iance under trans- 
lation along [111] in the reciprocal lattice is expressed as 

-^{^x ~^ "^'^2/ ^ "^^y^'^i ^ '^^^) ~ i^x y y ^z) ' This 

implies C{k'^) = C{k'^ + \/37r/3), or that the period of 
Chern number is one third of the period along k'^ . On the 
planes k'^ = and k'^ = \/37r/2, which are invariant under 
Mz : z^ —z\ one recovers the Cq symmetry, and we 
use Eq. (|3j) to evaluate their Chern numbers. Using the 
Ce eigenvalues at high symmetry points on k'^ = V^7r/2 
and k'^ = 0, we conclude that C{k'^ =0) = 6n — 2 and 
C{k^^ = A/37r/2) = 6n, which implies again the existence 
of Weyl bulk gapless points in the material, related by 
Cq * Mz' and lying along TZ' direction (see Fig[2fc,d)). 
At any /c^ the symmetry is {Cq * Mz)'^ = C3, which does 
not distinguish between C = —2 and (7 = 1 Chern num- 
bers, and as the difference in Chern numbers between 
k'^ = 0, V^7r/2 is —2, we could, in principle, have several 
gapless scenarios for /c^ > 0, the simplest of which are: 
1) a double- Weyl fermion at some k'^ = kc^k'^ = ky = 0] 
unlike the case of C4, this is not protected by symme- 
try, and will generically split into multiple linear Weyl 
fermions. 2) two linear Weyl fermions at two different 
/c^ > 0, each of which would give a change in Chern 
number —1; this is not possible in the current situa- 
tion due to the C3 eigenvalues of the two bands cross- 
ing: the Cs matrix is diagonal with e*^/^, e*^^/^ and 
can differentiate between C = —1 and C = — 2, disal- 
lowing the presence of the former. 3) four linear Weyl 
fermions, one of them with monopole charge +1 located 
at /c^ = /c^ = 0, /c^ = /Cci, and some energy Ei and three 
Weyl fermions of monopole charge —1, at the same en- 
ergy £^2, at momenta away from the [111] direction and 
related by C3 symmetry; this is the case present in the 
realistic material, see Fig [2]. The k • p Hamiltonian de- 
scribing the two bands which give rise to the four Weyl 
points is that of Eq|4] (double- Weyl point) augmented 
by a linear term /c+(7+ + /c_cr_. Depending on whether 
El > E2 or El < E2 (material parameter dependent), 
the Chern number progression starting from /c^ = to 
/c^ = TT is —2, —3, or —2, 1, as one successively crosses 
the Weyl points. Unlike in the [001] case, as Mz' is not a 
symmetry of the system, no stable nodal line is present. 

In conclusion, we have presented a theory relating the 
Chern number of 2D insulators with Cn rotational sym- 
metry to the eigenvalues of the PCS rotation matrix at 
high-symmetry points in the BZ. We then analyzed 3D 
materials with a Cn rotation axis and show that if the 
Chern number as a function of the momentum parallel 
to the rotation axis has a jump of 2, the 3D material 
can, in special situations, exhibit a double-vortex gap- 
less point protected (by the PGS) from splitting into two 
usual Weyl points. We use this theory to show that fer- 
romagnetic HgCr2Se4 in 3D possesses a quadratic cross- 
ing (in the kx^ ky plane) if the ferromagnetism is aligned 



5 




FIG. 3: Momentum resolved spectral weight (oc 
lm[G{ky, kz, E)]) calculated on the surface after the system is 
cut along the /cx- direct ion. The energy is varied as a function 
of kz to ensure that it is inside the bulk gap. The intensity 
is normalized such that unity means the state is completely 
localized on the surface. The inset shows the band crossing 
on rZ, where each band, from the second to the fifth, is 
labeled by its C4 eigenvalue ^. 

along the [001] direction, with the double- Weyl node oc- 
curring along the TZ line. In this case, the system also 
enjoys a mirror symmetry, whose change in eigenvalues 
along some line in the kx^ ky plane renders this line nodal. 
The electrical polarization inside this line is l/2eL per 
kx^ky point, while it vanishes outside, creating a topo- 
logical multiferroic. If the ferromagnetism in HgCr2Se4 
lies along the [111] direction, the system is invariant un- 
der a combined symmetry operation C6 * M^' with the 
characteristic double- Weyl node along the TZ' line un- 
protected and split into four linear Weyl points. 
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